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Abstract 



In this paper, we analyze the subspace migration that occurs in limited-view inverse scattering problems. 
Based on the structure of singular vectors associated with the nonzero singular values of the multi-static 
response matrix, we establish a relationship between the subspace migration imaging function and Bessel 
functions of integer order of the first kind. The revealed structure and numerical examples demonstrate 
why subspace migration is applicable for the imaging of small scatterers in limited- view inverse scattering 
problems. 
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1. Introduction 

This work focuses on subspace migration for the imaging of small perfectly conducting cracks located in 
the homogeneous, two-dimensional space R 2 in limited- view inverse scattering problems. Subspace migration 
is a simple and effective imaging technique, and is thus applied to many inverse scattering problems. For 

£T) ■ this reason, various remarkable pro perties have been investigated in many studies. Related work can be 

(N ; found in 0, [2, 0, B E3, El EH IH and the references therein. 



Based on these studies, it turns out that subspace migration can be applied not only in full-view but also 
limited-view inverse scattering problems. However, the use of subspace migration-based imaging techniques 
has been applied heuristically. Although the structure of subspace migration in full- view inverse scattering 
problems has been established [l|, [7|, it is still applied heuristically in limited- view problems. Hence, a 
mathematical identification of its structure in such problems still needs to be performed, which is the 
motivation for our work. 

In this paper, we extend the analysis j7[ of the structure of subspace migration in the full-view problem 
to the limited-view problem of imaging perfectly conducting cracks of small length. This is based on the 
fact that the far-field pattern can be represented by the asymptotic expansion formula in the presence of 
such cracks. From the identified structure, we find a condition for successful imaging. Subsequently, the 
structure of multi-frequency subspace migration is analyzed. This shows how the application of multiple 
frequencies enhances its applicability. 

This paper is organized as follows. In Section [21 we survey two-dimensional direct scattering problems, 
the asymptotic expansion formula in the presence of small cracks, and subspace migration. In Section [31 
we investigate the structure of single- and multi-frequency subspace migration in limited-view problems by 
finding a relationship with integer-order Bessel functions of the first kind. Section [5] presents some numerical 
examples to support our investigation. 
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2. Preliminaries 

In this section, we briefly introduce two-dimensional direct scattering problems in the presence of small, 
linear perfectly conducting crack(s), the asymptotic expansion formula, and subspace migration. 

2.1. Direct scattering problems and the asymptotic expansion formula 

Let E be a straight line segment that describes a perfectly conducting crack of small length 2£, 

E(z) := {z = {x,y) T : -£ < x < £} . 

With this E, let Utot(x, 8; k) satisfy the following Helmholtz equation 

J Au tot (x,0;fc) + fc 2 u tot (x,0;fc)=O in M 2 \E 
[ u to t(x, 9;k) =0 on E, 

where k = 2ir/\ denotes a strictly positive wavenumber with wavelength A and we assume that k 2 is not an 
eigenvalue of ([1]). 

It is well known that the total field u tot (x, 9; k) can be decomposed as 

u tot (x, 9; k) = w inc (x, 9; k) + M sca t(x, 9; k), 

where Ui nc (x, 9; k) = e lke ' x is the incident field with direction 9 n on the two-dimensional unit circle S 1 , and 
Mscatfx, 9; k) is the unknown scattered field satisfying the Sommerfeld radiation condition 

lim |x|Va ( ^scat(xe;fc) _ ^ \ = 

uniformly in all directions x = x/|x|. 

The far-field pattern Uoo(x, 0; k) of the scattered field u sca t(x, 9; k) is defined on S 1 and can be represented 

as 

e *fe|x| / /i 

u scat (x,0;fc) = ——j I Uoo(±,d;k) +0 ( i^T 

uniformly in all directions x = x/|x| and |x| — > +oo. Based on |J|, the far-field pattern can then be 
represented as the following asymptotic expansion formula. 

Lemma 2.1 (Asymptotic expansion formula in the presence of a small crack). Let w, to t(x, 9; k) satisfy flj) 
and Uinc(x, 9; k) = e l ' ce "' x . Then, for < I < 2, the following asymptotic expansion holds: 

2tt 



Uoo (x,9;k) = -— — u inc (z,6>;fc)M inc (z,x;fc) + 0( — -— J. (2) 

2.2. Introduction to subspace migration 

At this point, we apply @ to explain an imaging technique known as subspace migration. From [2j, 
subspace migration is based on the structure of singular vectors of the collected Multi-Static Response 
(MSR) matrix 



K(k) = 



^oo\X m , t/ n , Kj 



N 



m,n=l 



whose elements WoofXm, 9 n ; k) make up the far-field pattern with incident number n and observation number 
to. From now on, we assume that there exist S different small cracks E s with the same length £, centered 
at z s , s = 1, 2, • • • , S. The asymptotic expansion formula © can then be represented as follows: 



27T A , _ ., ; : rv 2tt 



Uoo (x,0;/c)«-— — ^ Uinc (z s ,6»„;fc)u inc (z s ,x m ;fc) = -— — ^ e ^( e " *™^ . (3) 



2.3. Introduction to subspace migration 

The subspace migration imaging algorithm is based on the structure of the singular vectors of the MSR 
matrix K(k). For the sake of simplicity, suppose that the incident and observation directions coincide. In 
this case, for each x m = — 9 m , the mn—th element of the MSR matrix becomes 

2vr 



(-B 9 -b\- y r »fc(e„+e m )-z a 

,{ 6 m ,V n ,k)- ln(£/2) Z. e 



Now, let us define a unit vector 



W(x; k) := — 






such that K(fc) can be decomposed as 

K(fe): 



2?r 



\n(£/2) ^ 



^W(z s ;fc)W(z s ;fc) 5 



Using this decomposition, we can introduce the subspace migration imaging algorithm. By performing 
Singular Value Decomposition (SVD) on K(k) 



K(Jfc) = U(fc)B(fc)V(fc)* w^CT s (fc)U s (fc)V s (fc)*, 



we can observe that the left- and right-singular vectors, {U s (fc)} s=1 and {V s (fc)} s=1 , respectively, satisfy 

U s (fc) ~ W(z s ;fc) and V s (k) ~ W(z s ; jfc). 
Hence, by defining an imaging functional 



J(x; k) := 



J2 (W(x; k), Ua(fc)) (W(x; fc), V.(fc)) 



we can obtain the shape of E s based on the orthogonality of the singular vectors. Here, (a, b) = a • b. A 
more detailed discussion can be found in [2J, [D, |9|, |ljj, lll|, [l2|, [l3| . 



3. Structure and properties of subspace migrations in limited-view problems 

In this section, we carefully explore the structure of subspace migration in limited view problems and 
discuss its properties. Suppose that each 6 n is an element of Sg ub C S 1 such that 

6 n := (cos6>„, sin0„) , 8„ = a + ((3 — a)- 



N -1' 

where a^O and j3 ^ 2n. From this, we will derive a useful approximation in the following Lemma. This 
will play a key role in our exploration. 

Lemma 3.1. Assume that k is sufficiently large, £ S 1 ^, and x, z £l 2 . Then, the following approximation 
holds: 



f e lke -^~ 2) dS{6) = {p- a)J (fc|x - z|) + 0[ —=L 



\/k\K-A 



Proof. Let us consider polar coordinates, such that = (cos 9, sin 9) and x z = r(cos0, sm0). Because 
the Jacobi- Anger expansion 

oo 

e lzcose = J (z) + 2^2i n J n {z)cos(n9), 

n=l 

holds uniformly (see [5(), we can evaluate 

rP rP r0 °° 

e lk9(x - z) dS(e)= e jfercos(e -* ) dS'(6>) = / J (kr)d9 + 2 ^i n J n {kr)cos{n9)d9 

Ja Ja Ja n=1 

oo „,3 

= (/3 - a)Jo(kr) + 2 ^ i n J n {kr) I cos(n9)d9 

n=l J° 

°° i n 1 / \ 1 / 

= (/? - a)Jo(kr) +4^ — J n {kr) sin- ( n(/3 - a) J cos - I n(/3 + < 



Assume that fc|x — z| is sufficiently large. As the following asymptotic form holds 

J„(fc|x - z|) = 
we can neglect the term 



7rfc|x — z| 



<<>s J k\x- z| — ) 1 =() 



\/k\x-z\ 



(1) 



4^ — J n (kr) sin- (n((3 -a) j cos-f n(/3 + a) 

n— 1 ^ ' ^ 



This completes the proof. 

Based on Lemma |3. 11 we can immediately obtain the following results. 
Theorem 3.2 (Single-frequency subspace migration). For sufficiently large N, X(x; k) can be written as 



a 



Z(x; fc) « V ( J (fc|x - z s |) + of -— L 
^V W fc l x 



(5) 



Proof. Since the incident and observation direction configurations are the same, we set A6 p := \8 P — #p-ij 
for p = 2, 3, • • • , AT, and A<?i := \6i — 8^\. Then, applying Lemma IBTTl yields 



X(x; k) 



J2 (W(x; k), U s (fc)> <W(x; fc), V s (fc)) 



S / N 

E E< 



: ikO p -{-x.— z s ) 



A8p 

2ir 



1 



^ik6-(x— z s ) 



(/3 ~ a) 2 W sl 

\ sul 

Therefore, we can obtain 



dS(0) 



03 



i^g^-aJJbWx-^J + O^ 



Z(x; fe) « £ f J (fc|x - z s |) + o(-—L_ 



This finishes the proof. □ 

Based on recent work (see [l|, y, |g, |3, [9|, [l0|, [llj, [l2j, [l3|), it has been confirmed that a multi- frequency 
approach gives better results than applying a single frequency. The following theorem supports this fact. 

4 



Theorem 3.3 (Multi-frequency subspace migration). Let kf := 2n/Xf and N be sufficiently large. Then, 
multi-frequency subspace migration 



X MF (x,fci,fcjr;.F) = 



1 



F S 



E E < w ( x - fc /)> u *( fc /)> < w ( x - fc /)> v ^/)> 

/=ls=l 



can be written as follows 



Xmf(x, ki,k F ;F) 



kf 



kp — k\ 



J (A;f|x — z|) 2 + Ji(/cf|x - z| 

fci 



fcj? — fci 



J (fci|x-z|) 2 + Ji(fci|x-z| 



(6) 



Proof. Similar to the proof of Lemma [3.11 we consider polar coordinates such that 6 — (cos 9, sin 8) and 
x — z s = r s (cos(j> s , sm(j) s ). Then, according to ([5]), we can observe that 



F S 



S ,fc, 



'dfe. 



Imp(x; fc) « i EE J «( fc l x ~ z ^l) 2 w fe-TT- E / ' J °( fcr *) 2 

/=ls=l F 1 s=l ^ fel 

Using this, we apply an indefinite integral formula of the Bessel function (see [1J, page 35] 

J Q (x) 2 dx = x( J (a;) 2 + Ji(a;) 2 J + / Ji(a;) 2 da;. 



Thus, we immediately obtain 

J (kr s ) 2 dk = k F (j (k F r s ) 2 + Ji(k F r s )A - k x fj (fcir s ) 2 + Ji(hr s )A + / ? Jx{kr 8 fdk. (7) 



Based on [7, Theorem 3.4], the last term of Q can be disregarded. Hence, we obtain ©. 



D 



Based on Lemma 3.1, Theorems 3.2 and 3.3, we can examine certain properties. They can be summarized 
as follows. 

1. The condition of sufficiently large k is a very strong assumption. If this condition is not fulfilled, the 
term 

01 ' 



V fc l x ^ z l 

of ([5]) will disturb the imaging performance. This is why the application of a high frequency yields 
good results, and why subspace migration in the full-view problem offers better results than in the 
limited- view case. 
2. If x is close to z, i.e., |x — z| is small, then fc|x — z| is not sufficiently large and we cannot use the 
asymptotic form (0]). This means that the remaining term of Lemma I3~T1 



°° i n 1 / \ 1 / 

4E — J n (kr s ) sin- (n([3 ~ a) I cos - I n((3 + a) 



is not negligible, so some blurring effects will appear in the neighborhood of cracks (refer to Figured]). 
3. Conversely, if the location of x is far from z, i.e., |x — z| is large, then the map of I(x; k) will be plotted 
at 0. 



4. Numerical examples 

In order to support some identified properties in Theorems 13.21 and 13.31 we now present some numerical 
examples. For this purpose, N = 12 different incident and observation directions are applied, such that 

77 — 1 7T )7T 

d n = a + (/? - a) _ for a = — and /3 = —. 

The applied wavenumber is of the form kf = 2ir/\f, where \f, f = 1, 2, • • • , F(= 10), is the given wavelength 
and kf are equi-distributed in the interval [ki, fcio] with Ai = 0.6 and Aio = 0.2. 

The far-field pattern data u^, (x; fc) are generated by a Fredholm integral equation of the second kind 
along the crack introduced in [8|, Chapter 4] to avoid committing inverse crimes. Three cracks E s with small 
length £ = 0.05 are chosen for the numerical simulations, with 

Ei = {[* - 0.6, -0.2] T : -t < t < £} 

E 2 = {£**/*[* + 0-4, t + 0.35] T : -£ < t < 1} 

S 3 = {5K W6 [t + 0.25,t-0.6] T : -£<t<£} . 

Here, 9t v denotes rotation by <p. 

Figure [1] shows maps of X(x, fcio) and Imf(x, k\, fcio;10) for the full- and limited-view cases. Based on 
Theorem 13.21 the locations of the E s are successfully identified, but many replicas and blurring effects in 
the neighborhood of cracks disturb the identification. However, applying a multi-frequency yields a more 
accurate result than the single-frequency case. This supports Theorem l3.3l In any case, we can conclude that 
our analysis demonstrates that subspace migration can determine the location of small cracks in limited- view 
problems. 

5. Conclusions 

In this paper, we considered subspace migration for the imaging of small, perfectly conducting cracks. 
Based on a relation with an integer-order Bessel function of the first kind, we confirmed the effectiveness of 
subspace migration in limited- view inverse scattering problems under certain strong assumptions. We also 
identified some particular properties of the structure of this migration. 

The main subject of this paper is the imaging of small cracks. Extension to the imaging of arbitrary 
shaped, arc-like cracks will be considered in future work. Although we considered a two-dimensional problem, 
its extension to three dimensions will be an interesting challenge. 

Acknowledgements 

This research was supported by the Basic Science Research Program through the National Research 
Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (No. 2012- 
0003207) , and the research program of Kookmin University in Korea. 

References 

[1] H. Ammari, Mathematical Modeling in Biomedical Imaging II: Optical, Ultrasound, and Opto-Acoustic Tomographies, 

Lecture Notes in Mathematics, 2035, Springer- Verlag, Berlin, 2011. 
[2] H. Ammari, J. Gamier, H. Kang, W.-K. Park, and K. S0ma, Imaging schemes for perfectly conducting cracks, SIAM J. 

Appl. Math. 71 (2011) 68-91. 
[3] H. Ammari, E. Iakovleva and D. Lesselier, A MUSIC algorithm for locating small inclusions buried in a half-space from 

the scattering amplitude at a fixed frequency, SIAM Multiscale Model. Sim. 3 (2005) 597-628. 
[4] H. Ammari, H. Kang, H. Lee and W.-K. Park, Asymptotic imaging of perfectly conducting cracks, SIAM J. Sci. Comput. 

32 (2010) 894-922. 
[5] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Springer, New York, 1988. 

6 




ro 




(a) Map of Z(x,0.2) 



(b) Map ofZ(x,0.2) 




u 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




-1 -0.5 0.5 1 

x-axis 

(c) Map of 1mf(x, fci, fcio; 10) (d) Map of Imf(x, fei, fcio; 10) 

Figure 1: Maps ofZ(x, fcio) and 1mf(x, k\, fcio; 10). Left column: full-view case. Right column: limited- 



view case. 



[6] S. Hou, K. Huang, K. S0lna, and H. Zhao, A phase and space coherent direct imaging method. J. Acoust. Soc. Am. 125 

(2009) 227-238. 
[7] Y.-D. Joh, Y. M. Kwon, J. Y. Huh and W.-K. Park, Structure analysis of single- and multi-frequency subspace migrations 

in inverse scattering problems, Prog. Electromagn. Res. 136 (2013) 607-622. 
[8] Z. T. Nazarchuk Singular Integral Equations in Diffraction Theory, Karpenko Physicomechanical Institute, Ukrainian 

Academy of Sciences, Lviv, 1994. 
[9] W.-K. Park, Analysis of a multi-frequency electromagnetic imaging functional for thin inclusions, in reivision, available at 
http : //arxiv . org/abs/1208 . 2063 
[10] W.-K. Park, Non-iterative imaging of thin electromagnetic inclusions from multi-frequency response matrix, Prog. Elec- 
tromagn. Res. 106 (2010) 225-241. 
[11] W.-K. Park, On the imaging of thin dielectric inclusions buried within a half-space, Inverse Problems 26 (2010), 074008. 
[12] W.-K. Park and D. Lesselier, Fast electromagnetic imaging of thin inclusions in half-space affected by random scatterers, 

Waves Random Complex Media 22 (2010) 3-23. 
[13] W.-K. Park and T. Park, Multi-frequency based direct location search of small electromagnetic inhomogeneities embedded 

in two-layered medium, Comput. Phys. Commun. 184 (2013) 1649-1659. 
[14] W. Rosenheinrich, Tables of Some Indefinite Integrals of Bessel Functions, availabe at 

http : //www . f h- j ena . de/~ rsh/Forschung/Stoer/besint . pdf 



